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Abstract: A simple ansatz is suggested for the structure of threshold resummation of 
the momentum space physical evolution kernels ('physical anomalous dimensions') at all 
orders in (1 — x), taking as examples Deep Inelastic Scattering and the Drell-Yan process. 
Each term in the expansion is associated to a distinct renormalization group and scheme 
invariant perturbative object ('physical Sudakov anomalous dimension') depending on a 
single momentum scale variable. Both logarithmically enhanced terms and constant terms 
are captured by the ansatz at any order in the expansion. The ansatz is motivated by 
a large-/3o dispersive calculation. A dispersive representation at finite (3q of the physical 
Sudakov anomalous dimensions is also obtained, associated to a set of 'Sudakov effective 
charges' which encapsulate the non-Abelian nature of the interaction. It is found that the 
dispersive representation requires a non-trivial, and process-dependent, choice of variables 
in the (x, Q 2 ) plane. Some interesting properties of the physical Sudakov anomalous dimen- 
sions are pointed out. The ensuing 1/N expansion in moment space is straightforwardly 
derived from the momentum space expansion. 
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1. Introduction 

Threshold resummation, namely the resummation to all orders of perturbation theory of 
the large logarithmic corrections which arise from the incomplete cancellation of soft and 
collinear gluons at the edge of phase space, is by now a well developed topic [1, 2] in 
perturbative QCD. About ten years ago, the subject was extended [3-5] to cover also the 
resummation of logarithmically enhanced terms which are suppressed by some power of 
(1 — x) for x — ► 1 in momentum space (or by some power of 1/N, N — > oo in moment 
space), concentrating on the case of the longitudinal structure function Fl(x, Q 2 ) in Deep 
Inelastic Scattering (DIS) where these corrections are actually the leading terms. As far 
as I am aware, little work has been performed on this subject since then. In this paper, 
building upon the recent work [6], I provide a very simple ansatz for the structure of 
threshold resummation at all orders in (1 — x), working at the level of the momentum 
space physical evolution kernels (or 'physical anomalous dimensions', see e.g. [7-10]), which 
are infrared and collinear safe quantities describing the physical scaling violation, where 
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the structure of the ansatz appears to be particularly transparent. I shall deal explicitly 
with the examples of the (non-singlet) Deep Inelastic Scattering (DIS) structure function 
F2(x,Q 2 ), as well as with the Drell-Yan process. The ansatz is motivated by a large-/?o 
dispersive calculation, and, following [6,11-15], easily generalizes itself to a general finite- 
Po dispersive representation of the physical evolution kernel at any order in the x — > 1 
expansion. 

The paper is organized as follows: section 2 is devoted to DIS. In section 2.1, the 
momentum space ansatz for threshold resummation for the non-singlet structure function 
F2(x,Q 2 ) is displayed, which introduces at each order of the x — > 1 expansion a new 
'jet physical anomalous dimension'. The ansatz is justified in section 2.2 by a large—/?o 
calculation, which also provides a dispersive representation for each of the previous physical 
anomalous dimensions. These dispersive representations are then extended to finite-/3o in 
section 2.3, where a set of 'jet Sudakov effective charges' which encapsulate the non-Abelian 
nature of the interaction is introduced. The Drell-Yan case is addressed in section 3 in quite 
a similar way. The ansatz for the analogous r — > 1 expansion is given in section 3.1, and 
justified in section 3.2 by a large-/3o calculation, which yields the dispersive representation 
of the corresponding 'soft physical anomalous dimensions' specific to the Drell-Yan process. 
These representations are then extended to finite-/3o in section 3.3, where a set of 'soft 
Sudakov effective charges' is introduced. Both in the DIS and in the Drell-Yan case, a non- 
trivial (process-dependent) choice of expansion parameter in the (x,Q 2 ) (resp. (r, Q 2 )) 
plane has to be made in order to derive the dispersive representations. The conclusions are 
given in Section 4. The 1/N expansion in moment space is derived in a straighforward way 
from the corresponding momentum space expansion ansatz in Appendix A (DIS) and B 
(Drell-Yan). These two last appendices also clarify the connection between the definitions 
of the 'jet' (DIS) (or the 'soft' (DY)) scales in momentum and moment spaces respectively. 

2. Deep Inelastic Scattering case 
2.1 A systematic expansion for x — > 1 

The scale-dependence of the (flavour non-singlet) deep inelastic structure function F2 can 
be expressed in terms of F2 itself, yielding the following evolution equation (see e.g. Refs. [7, 



K(x,Q 2 ) is the momentum space physical evolution kernel, or physical anomalous dimen- 
sion; it is renormalization-group invariant. In [6], using known results of Sudakov resum- 
mation in moment space, the result for the leading contribution to this quantity in the 
x — > 1 limit was derived. For completness, I reproduce this short derivation here. 
Defining moments by 



9,10]): 




(2.1) 




(2.2) 
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Eq. (2.1) implies that the moment-space physical evolution kernel is: 



K(N,Q 2 ) = I dxx N - 1 K(x,Q 2 ) 



dlnF 2 (N, Q 2 ) 
dlnQ 2 



(2.3) 



Let us now consider the N — > oo limit (corresponding to x — > 1 in momentum space). 
In this limit the evolution of the structure function takes a simple from [11, 12, 16, 17] (a 
straightforward derivation from the standard Sudakov resummation formulas can be found 
in [14]): 



dlnF 2 (N,Q 2 



dx 



x 



N-l 



1 



J ((1 - x)Q 2 ) + H(a s (Q 2 )) + 0(1/N) (2.4) 



dlnQ 2 jo 

where the first term, which includes the N — > oo divergent corrections to all orders, is 
controlled by the 'jet' physical Sudakov anomalous dimension J([jl 2 ) (a renormalization 
scheme invariant quantity): 

dB(a s (n 2 )) 



J{^)=A{a s {^)) + 



(2.5) 



din ix 2 

where A(a s ) is the 'cusp' anomalous dimension, and B(a s ) the standard 'jet' Sudakov 
anomalous dimension 1 . Moreover, the constant term can be written [14] in terms of the 
quark electromagnetic form factor T(Q 2 ) [18-21]: 



H{a s {Q 2 )) 



din {HQ 2 ))' 
dlnQ 2 



+ 



f 

Jo 



dfj 2 
V 2 



(2.6) 



(l,a s (QV = 0) + B(a s (Q 2 )) , 



where each of the two terms in the first line is separately infrared divergent, but the 
divergence cancels [14] in the sum; in the second line the result is expressed in terms of 
B(a s ) and G (Q 2 /^ 2 , a s ( / u 2 ), e) , which is the finite part of din [F(Q 2 )) 2 /din Q 2 as defined 
in Ref. [21] using dimensional regularization. 

Comparing (2.3) and (2.4) one therefore finds the following relation [6] in momentum 
space: 



K(x,Q 2 ) = 



j((l-x)Q 2 ) | dln{F(Q 2 )f 
1 — x 



+ 



J ((1 - x)Q 2 ) 



+ 



+ 



dlnQ 2 

^ln(^(Q 2 )) 2 ( fdj^ J } 
Jo V 2 



dlnQ 2 



+ 



V 



6(1 -x) 



(2.7) 



infrared finite 



7 



+ 0((l-x) ) , 
where the integration prescription []+ is defined by 



dx F(x) 



J{(1-x)Q 2 ) 



1 - x 



j\( m - ni )(£«£m\, ( ,8, 



1 A and B are separately scheme dependent quantities. 
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where F(x) is a smooth test function. This prescription accounts for the divergent vir- 
tual corrections, which cancel against the singularity generated when integrating the real- 
emission contributions near x = 1. One thus finds that J ((1 — x)Q 2 } /(l — x) is the leading 
term in the expansion of the physical momentum space evolution kernel K(x,Q 2 ) in the 
x — > 1 limit with (1 — x)Q 2 fixed. The term proportional to 5(1 — x) is comprised of purely 
virtual corrections associated with the quark form factor. This term is infrared divergent, 
but as indicated in the second line in (2.7), the singularity cancels exactly upon integrating 
over x with the divergence of the integral of J ((1 — x)Q 2 ) /(l — x) near x — > 1. 

Next I observe that eq.(2.7) strongly suggests the following generalization to a sys- 
tematic expansion for x — > 1 in powers of 1 — x, or, more conveniently (for reasons to be 
clarified in section 2.2), in powers of 

r = 1— ^ (2.9) 

x 

at fixed jet mass 

W 2 = rQ 2 (2.10) 

namely: 

Q 2 ) = \j (W 2 ) + dln J^ 2 2))2 «(1 " x) + J (W 2 ) + r J l (W 2 ) + O (r 2 ) , (2.11) 

where (barring the virtual contribution) all coefficients J (W 2 ) and Ji (W 2 ) are renor- 
malization group and scheme invariant 'effective charges' [7], the physical 'jet' Sudakov 
anomalous dimensions, functions of a single variable-the jet mass W 2 , that can be com- 
puted order by order in a s (W 2 ). Eq.(2.11) represents a very simple momentum-space 
version of Sudakov resummation for the physical evolution kernel. I note that each power 
of r could be a priori multiplied by some powers of lnr. However, evidence from a large-/?o 
calculation (see below) suggests this is actually not the case, and only powers of r even- 
tually do appear. According to the ansatz (2.11), the whole towers of Sudakov logarithms 
ln p (r)/r in leading order, as well as those which are suppressed by a power of r (rMn J (r), 
i > 0) in subleading orders, are thus generated by expanding the J (W 2 ^ and Ji (W 2 } 
physical jet anomalous dimensions in powers of a fixed (i.e. x-independent) coupling, say 
a s (Q 2 ), i.e. are ultimately generated by standard renormalization group logarithms. In 
particular, any given term in this r expansion is associated to an infinite tower of arbi- 
trary powers of lnr. Moreover, the same ansatz captures all non-logarithmic terms, i.e. 
all In (r)/r terms in leading order (together with the virtual contribution), as well as all 
terms r % ln°(r) [i > 0) at subleading orders . The conjecture that all Sudakov logarithms in 
the physical evolution kernel can be absorbed by the change of argument Q 2 — > W 2 in the 
running coupling was already proposed in [7], a generalization of the suggestion in [22,23]. 
In the leading order of the r expansion, this statement is now seen to be a consequence 
of the standard resummation formalism [1,2]. Beyond leading order, this statement re- 
mains a conjecture, which could however be checked by matching the ansatz (2.11) with 
existing [24, 25] fixed order calculations; moreover, this comparison would determine the 
perturbative expansion of the J {W 2 } and Ji (W 2 ^ anomalous dimensions up to 0(a^). 
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I also note that no integration over the Landau pole explicitly appears in this momen- 
tum space version of Sudakov resummation. We shall see however that the perturbative 
expansions of the subleading physical Sudakov anomalous dimensions Ji (W 2 ) do contain 
infrared renormalons (at the difference of J (W 2 )l). 

One can also give the generalization of eq.(2.7), with a regularized virtual contribution: 

K(x,Q 2 ) = ^J{W 2 ) + ^^Q? 2 S(l-x) + J (W 2 ) +rJ 1 (W 2 ) + O (r 2 ) 



- J (W 2 ) 



( 



+ 



din {HQ 2 ))' 
dlnQ 2 



V 



f Q2 d^ 2 „, 2 
Jo A* 



\ 



5(1 -x) 



(2.12) 



l-x 



W 2 



1-1/12 



Q 



infrared finite 

+ Jo(W 2 ) +r Ji (W 2 ) +0(r 2 ) , 

where the integration prescription []+ is now defined by (r 

f 1 dxF(x) \- J (W 2 )} = ^ dx (F(x) - J (W 2 ) -F(l) - J— J ((1 - x)Q 2 ) ) , (2.13) 
Jo l r 1+ Jo ^ r t — X / 

where F(x) is a smooth test function. The proof that (2.12) is indeed equivalent to (2.7) 
for the two leading terms is given in Appendix A, by considering the moments of the two 
expressions. 

2.2 Dispersive representation of the physical Sudakov anomalous dimensions 
(large f3 ) 

The correctness of the ansatz Eq.(2.11) can be checked at large—/?o, upon taking the x — > 1 
expansion of the large-/?o dispersive representation of the physical evolution kernel. This 
procedure will actually yield a more powerful result, namely the (large-/?o) dispersive rep- 
resentations of the physical Sudakov anomalous dimensions J and J; L . 

Let us consider the large-/?o but arbitrary x (0 < x < 1) dispersive representation [8] 
of the physical evolution kernel. At the level of the partonic calculation, the convolution 
on the r.h.s of (2.1) becomes trivial in the large-/3o limit since the 0(a s ) corrections to 
F2(z,Q 2 ) generate terms that are subleading by powers of 1/Po- Therefore, in this limit 
K(x,Q 2 ) is directly proportional to the partonic dF 2 (x, Q 2 )/dln Q 2 , so 



K(x,Q 2 )\ 



large /3 



C F 



I 

Jo 



o A* 



dfi 2 
—ay 



"iihsfUriif.x). 



(2.14) 



where J-(/j, 2 /Q 2 ,x) is the standard notation for the characteristic function corresponding to 
F 2 (x, Q 2 )/x (see Eq. (4.27) in [8] or (3.2) in [26]) (with T = -p 



2 d 
dfj, 2 



and a^ ink (ji 2 ) is the 



integrated time-like discontinuity of the one-loop V-scheme coupling (which corresponds 
to the single dressed gluon propagator, using 'naive non-abelization'): 

V) . „ Min w ^ _ r dm 2 

dln^' ay (M ) = ~L ^ 



J Mink, 



-p v (m 2 ), 



(2.15) 
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where 

p v (n 2 ) = ^Im{ a V(-fi 2 -iO)/Tr} (2.16) 



and 



(2.17) 



vr fa In {k 2 /A 2 v ) ' 
with 

A^ = A 2 e 5/3 , (2.18) 

where A 2 is defined in the MS scheme, and (3q = j^Ca — gAfy. 

Next, one takes the x — ► 1 expansion under the integral (2.14) with a fixed invariant 
jet mass W 2 = Q 2 {1 — x)/x = Q 2 r. To achieve this, one splits the characteristic function 
into its real and virtual contributions: 

F(e,x) =J r{r) {e,x)e(l-x-ex)+Vs(e)5{l-x) (2.19) 

where e = fi 2 /Q 2 , and expand the real contribution (the virtual contribution should be left 
unexpanded). Using the explicit expression 2 for jrt r )(e,x) in [8], one obtains the small r 
expansion at fixed £ : 

xJ* r \e,x) = -^)(0 +^ r) (0 + r^ r) (0 + 0(r 2 ) , (2.20) 



where 



with 



4 r \0 = ln£ + I + 2£ln£ - 4£ + \ £ 2 (2.22) 

^ (r) (e) = -^ine-6-9£ine + ie-£ 2 ine + ^e 2 . 



Thus (• ee -y 2 *): 



with 



x& r \e,x) = -^ r) (0+H r) (0 + r^ r) (0 + O(r 2 ) , (2.23) 



■^ r) (0 = -l-2£ln£ + 2£-£ 2 (2.24) 
■^i (r) (O = ^ + 9ein£+ye + 2£ 2 ln£-10e 2 , 



2 The present normalization of J-(e,x) is half the one in [8]. 
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and, taking a second derivative: 



with 



fi r) (0 = -9£ln£ - ^ £ " 4^ 2 InC + 18 £ 2 
Using these results in (2.19), and noting that 

0(1 -x- ex) = #(£< 1) , 

one gets: 



35 



xF(e,x)= [-^7(0+^7o(0+»-^7i(0+0(r 2 )j + V s (e) 5(1 - x) , 



with 



^7(0=^ r) (0^<i) 
^(0=^ r) (Ofl«<i), 



(2.25) 



(2.26) 



(2.27) 



(2.28) 



(2.29) 



with 



xF{e,x)= [-^7(0+^7o(0 + »-^7i(0 + 0(r 2 )j + V s (e)<5(l-x) , 



^7(0=^ r) (0^<l) 
j>,(0=jf(0*K<l), 



(2.30) 



(2.31) 



and 



rl 



xF{e,x)= -jr j{0+ jr Jo{0+r jr Jii0+0{r ^ + Va (e) 5(1 - *) , 



(2.32) 



with 



^7(0=^ r) (0^<i) 
^(0 = ^ (r) (0^<i) , 



(2.33) 



where I used the fact that all the terms F^ft), and ^ r >(0, in (2.20) and 

(2.23) vanish at £ = 1 (which allows to treat the 9 function effectively as a multiplicative 
constant when taking the derivatives). These features actually follow from the stronger 
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property that the exact function J-( r \e,x), as well as its first derivative j-( r \e,x), both 
vanish at £ = 1, i.e. for e = Indeed one finds, expanding at fixed r: 

2 ;^)(e,x) = ^-(l + r)(3 + r)(l-0 2 + O((l-0 3 ) • (2-34) 

I further note that the terms which vanish for £ — > in ^j. (but not in TjX) are 
arithmically enhanced, hence non-analytic, which implies (see below) that the subleading 
physical Sudakov anomalous dimensions Ji do have renormalons (at the difference of the 
leading physical anomalous dimension J\). 

Reporting the result Eq. (2.32) into (2.14), one thus obtains the small r expansion 
(r = (1 — x)/x) of the physical evolution kernel: 



+ r r^ a ^V)^(6 + <V)} (2-35) 
Jo A* J 

+ <5(l-x)C F f° ^a{^V)V s (e) ■ 
Jo A 1 



/o A 4 ' 

Comparing with (2.11) I deduce (£ = /U 2 /VF 2 ): 



.2 

/o A* 
and also (e = / u 2 /Q 2 ): 

,2 



(2.36) 



din (^(Q 2 )) 2 



dlnQ 2 



large /3 



= C> r ^£#*"V)V.(e) , (2.37) 
jo A* 



which checks Eq.(2.11) in the large-/?o limit, and in addition gives the (large-/?o) dispersive 
representations of the physical Sudakov anomalous dimensions J and J, as well as the 
(large— /3o) dispersive representation of the quark form factor [14]. I note that the integral 
on the r.h.s. of (2.37) is ultraviolet convergent, but infrared divergent, as expected from the 
general properties [18-21] of the quark form factor. Indeed from the explicit expression in [8] 
one gets: V s (e) = C(lne/e) for e — > 00, whereas V s (e) — ► —1 for e — ► 0. For completness, I 
also give the dispersive representation of the regularized virtual contribution of Eq. (2.12). 
After a little algebra one finds, using (2.36) and (2.37): 



rf '° ( f ( g ))2 + f d 4j<s) = c F r d 4av>-v) (v.(,vq 2 ) + wiq 

« m V jo A 4 jo A 4 v 

where the integral on the r.h.s. is indeed both infrared and ultraviolet convergent. 



2 ) 

(2.38) 
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2.3 Dispersive representation of the physical Sudakov anomalous dimensions 
(finite (5q) 

Following [6,13], it is straightforward to give the generalization of Eqs. (2.36) and (2.37) 
at finite (5q: 



Jo M 9 (2.39) 
r°° fin 2 

^(W 2 ) = C F % a ^r\^)Pj^), 
Jo M 

where, in full analogy with (2.15), 

^)-^ : ^>-/;t^>> , 2 .4o, 

and, similarly to (2.16), pj,j x {^ 2 ) correspond to the timelike discontinuities of some "Eu- 
clidean" effective charges, originally defined for spacelike momenta: 

Eucl i 1 2 \ _ f d/J, Mink , 2 \ fJ> /k f°° dfT 2 

*J,jfr)=] Q — ajiJi (j.) il + ii2/k2)2 =-J Q ^PW,^)- (2-41) 



Furthermore, the generalization of (2.37) is (e = fi /Q 



2/n2v 



dln(HQ 2 )) 2 _ n rdn 2 nMinkf ^.. : 



d In Q 2 J /J, 2 

The generalization of (2.38) is: 



= C F ^-af m \^)V s {e) . (2.42) 



+ af nk { i a 2 ) (Pji^/Q 2 )-!) 



(2.43) 



where terms have been arranged properly to have both infrared and ultraviolet conver- 
gence: I note that the integrals over cijj mk ([i 2 ) and aj mk (fL 2 ) in (2.43), although infrared 
convergent, are separately ultraviolet divergent. 

From these dispersive representations, the perturbative expansions of the 'jet' Sudakov 
effective charges aj mk (/j, 2 ) and aj. mk (fi 2 ) can be derived in a straightforward way [6] 
order by order in the full non-Abelian theory, given the expansions of J(W 2 ) and Ji{W 2 ) 
(and similarly for a^ ink (fi 2 ) using (2.43) and (2.6)). 

3. Drell-Yan case 

3.1 The r — ► 1 expansion 

The cross section of the Drell-Yan process, h a + — > e + e~ + X, is: 
da _ 4:ira 2 m sr^ f 1 dxi dxj 



dQ 2 9Q 2 S J JQ .1 ; -' i 



^ r ^ doc ' dx ' 

22 / — — fi/h a (Xi,fJ<F) fj/h b (Xj,HF)gij{T,Q 2 ,H 2 F ) (3.1) 
. . Jn Xi Xj 
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where s is the hadronic center-of-mass energy, s = XiXj-s is the partonic one, Q 2 is the 
squared mass of the lepton pair and r = Q 2 /§. The partonic threshold r — ► 1 is character- 
ized by Sudakov logarithms. 

Let us define the Mellin transform of the quark-antiquark partonic cross section in 
(3.1), e.g. in electromagnetic annihilation g qq (r,Q 2 , fx F ) = e 2 [5(1 — r) + 0(a s )], by 

G qq (N, Q 2 ,n 2 F ) = f 1 dr t n ~ 1 g qq (r, Q 2 ,fi F ) . (3.2) 
J o 

In the TV — > oo limit one can derive (see e.g. [14]) the analogue of Eq. (2.4): 



d\nG qq (N,Q 2 ,v? F ) _ 



dlnQ 2 



drr JV - 1 K DY (r, Q z ) 
(3.3) 
25 (Q 2 (l - r) 2 ) + H DY (o> s (Q 2 )) + 0(l/iV) , 



1 r 7V-l _ 1 

dr— 

1 — r 



where K dy (t, Q 2 ) is the momentum space physical Drell-Yan evolution kernel defined 
for arbitrary r, and the iV-dependent terms are controlled by the 'soft' Sudakov physical 
anomalous dimension S: 



S(n 2 ) = A{a s {^)) + - 



2 „ , ldV{a s (n 2 )) 



(3.4) 



2 dln/i 2 

where V is the standard (scheme-dependent) 'soft' Sudakov anomalous dimension relevant 
to the Drell-Yan process. Moreover, the constant term can be expressed in terms of the 
analytically-continued electromagnetic quark form factor [14]: 

d\n\(F{-Q 2 ))\ 2 , ( Q2 dfj 2 
o V 2 



H DY ( as (Q 2 )) = ^ IV ;^ " X + / ^V), 



dlnQ 2 

where the infrared singularities cancel [14] in the sum, as in (2.6): 

din | (H-Q 2 )) | 2 , d M 2 2 

H UY (a s {Q )) = L TfT7^ + / — ) 



(3.5) 



dlnQ 2 7 /i 2 

(l,« s (Q 2 ),e = 0) +/3 (a s )^(a s (Q 2 ))l +ip(a s (Q 2 )) 



(3.6) 



with ft(a s (Q 2 )) = In 



. Thus, in momentum space we have [6]: 
2\\ | 2 

' ) ) + » m IW 

1 - r 
5(Q 2 (l-r) 2 )' 

+ 



, , N 5(Q 2 (l-r) 2 ) din | (^(-Q 2 )) | 2 / „ 
K DY (r, Q 2 ) = 2 l V , ^ ' + '^ Q 7 ^ 5(1 - r) + ((1 - r)°) 



1 -r 



(3.7) 



+ 



dlnQ 2 



o p 



-S m {n 2 ) 



5(l-r)+0((l-r) ) 



infrared finite 



7 
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We see that the physical Sudakov anomalous dimension S controls the leading term in 
the expansion of the momentum space physical Drell-Yan kernel (3.3) near threshold. The 
t — ► 1 limit is taken such that E DY = Q(l — r), corresponding to the total energy carried by 
soft gluons to the final state, is kept fixed. The virtual contribution, proportional to 6(1— r), 
is determined by the quark form factor, analytically-continued to the time-like axis. This 
term is infrared singular, but upon performing an integral over r this singularity cancels 
with the one generated by integrating the real-emission term 25 (Q 2 (l — r) 2 ) /(l — r) near 

T -► 1. 

Similarly to the DIS case, Eq. (3.7) suggests a generalization to an expansion for r — ► 1 
in powers of 1 — r at fixed E DY . As in the DIS case, however, it turns out (see section 3.2) 
that in order to derive a dispersive representation such an expansion is not the appropriate 
one. Instead, one should consider an expansion in powers of the alternative variable: 



1 - r 



(properly normalized so that r E 



1 — r for r — > 1) at fixed W 2 Y , with: 



(3.8) 



W 2 = r 2 Q 2 

r " DY — DY *V 1 



(3.9) 



namely: 



K dy (t,Q 2 ) = — s(w£ y ) + 



din | [H-Q 2 )) I 



dlnQ 



, S(1-t)+ S (W 2 y ) +r DY S l (W 2 Y ) +0 ( 

^*dy) 
(3.10) 

where the coefficients S and Si are the physical 'soft' anomalous dimensions appropriate 
to the Drell-Yan process. After regularizing the virtual contribution, Eq. (3.10) can be 
written as: 



( 



6(1 - r) + S (W 2 Y ) + r DY S 1 (W 2 Y ) + O (r DY ) 



' DY 



dln|(.F(-Q 2 ))|' 



dlnQ 2 



+ , d 4^ 

o A 1 



<5(l-r) 



infrared finite 

+ So (W 2 Y ) + r DY S 1 (W 2 Y ) + O (r 2 Y ) , 



(3.11) 



where the integration prescription []+ is defined (similarly to (2.13)) by (r DY = 2 



l-r 



W 2 Y = r 2 DY Q 2 ): 



f 1 dTF ( T ) \ J_ 5 1 = Z 1 |f(t) — 5 (W D 2 Y ) - F(l) -L- S ((1 - r) 2 Q 2 ) 

Jo L r DY J + Jo L r DY i — r 

(3.12) 

where F(t) is a smooth test function. 
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3.2 Dispersive representation of the physical Sudakov anomalous dimensions 
(large (3 ) 

In quite a similar way to the DIS case, we start from the large-/3o but arbitrary r (0 < r < 1) 
dispersive representation of the physical Drell-Yan momentum space evolution kernel, 

^ DY (r,Q 2 )| largc/3o =C F ^ OO ^py(^ 2 )(^ DY (/i 2 /Q 2 ,r)-^ DY (0,r)) , (3.13) 

where ^ r DY (/i 2 /(5 2 ) t) is the characteristic function corresponding to g q q(r, Q 2 , fif^). Here 
one has to use [6] this version of the dispersive representation involving the first derivative 
of the characteristic function, the second derivative being too singular, as we shall see, on 
the phase-space boundary. 

Let us first take the r — > 1 expansion under the integral while fixing the total energy 
radiated into the final state, £ DY = Q(l — r). Starting from: 

T DY {e,r) =4 r Y ) (e,r)0((l-r(l + v ^) 2 ) + V t (e) <5(1 - r) , (3.14) 

where e = fi 2 /Q 2 , and expanding the real contribution using the explicit expression 3 for 
f${e,r) in [8] (Eq. (4.87) there), one obtains: 



+ O(l-r). 

(3.15) 

I note that beginning at next-to-leading order (the C((l — r)°) contribution), the expansion 
becomes singular at 4/j, 2 /E 2 Y = 1, which corresponds to the phase-space boundary in the 
r — > 1 limit. Although at 0((1 — r)°) the singularity is still 4 integrable, it will become 
non-integrable after taking one derivative will respect to fi 2 and inserted into the integral 
on the r.h.s. of Eq. (3.13)! 

The solution to this difficulty consists in performing a change of variables. Namely, instead 
of 1 — r and fi 2 /E 2 Y , one should use r DY = 2 and 

' DY DY 

and perform the r DY — ► expansion in powers of r DY , at fixed £ DY . The result is: 

4 r Y } (e, r) = — (£ DY ) + ^2 Y (&v) + r DY ^g Y (£ DY ) + O(^dy) , (3-17) 



3 The present normalization of Tuy is half the one in [8]. 

4 Non-integrable singularities start appearing at order 0(1 — r). 
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with 

= 2tanh- 1 (Vl-4^ DY ) 
?0%(M = v/1-4Cdy-(1 " 4C DY )tanh- 1 ( v / l-4^ DY ) 

^3v(W = I ( ~ 19+ ^^ 2 — + ( 18 + 16 &~ + 32&) Umh-Hy/T^U)) 



(3.18) 



where all the terms T^(e,r) an d ^d Y Kdy) now vanish at the phase-space boundary 
4£dy = 1. I actually checked the stronger exact property that ^ Y (e,T) = for 4£ DY = 1, 
which explains the previous facts. Namely one finds, expanding at fixed r DY for 4£ DY — > 1: 

^(e,r) = ^ (l + J4) ^1 + ^dy V1-4^dy + ... • (3-19) 



-(r) 



Taking one derivative (' = — /x ^j) one thus gets the expansion: 



with 



^dy( £ i r ) — dy Kdy) + J 7 3 iDY (£dy) + TtsyF i jDY (£dy) + C(^dy) > 

''dy 



DY ) 



(3.20) 



•^0,dy(£dy) 
•^1,dy(^dy) 



=h= -^1-4^- 4^ DY tanh~ 1 ( A /l-4e DY ) 



9-38 g DY + 64 ^ 

8vi-4e DY 

r(r) t 



(3.21) 



2^(1 + 4e DY )tanh- 1 (Vl 3 4^ 



where all the terms ^dy((dy) an d ^ y (^dy) are now singular, but integrable, at the phase- 
space boundary 4£ DY = 1. This result follows immediately from (3.19) which yields: 

2 



1 + -7Y 2 



1 1 

2 VI - 4£dy 



Using these results in (3.14), and noting that: 

0((l-r(l + v ^) 2 ) = 0(4£ dy <1) , 

one gets: 



(3.22) 



(3.23) 



with 



'"dy 



■ Fs{M + ^So(^dy) + r DY F Sl (£dy) + 0(r 2 DY ) + V t (e) 5(1 - r) , (3.24) 



F Si (U) = ^2y(&y) ^ (4 ^DY < 1) , 



(3.25) 
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and 



^D Y (e,r) = 
with 



MM + M(M + r DY P Sl (M + 0{rl Y ) + V t (e) 5(1 - r) , (3.26) 



^s(^d Y ) = 4 T Y ) (eDY)^4e D Y < i) 

^,fey) =^%(U)0(±U <1) 



(3.27) 



where I used the previously mentioned fact that the terms on the r.h.s. of (3.17) all 
vanish at 4£ DY = 1 (which allows to treat the function effectively as a multiplicative 
constant when taking the derivatives). I further note that the terms which vanish for 
£dy — ► in Ts i (but not in J-^!) are logarithmically enhanced, hence non-analytic, which 
implies (see below) that the subleading physical Sudakov anomalous dimensions Si do have 
renormalons (at the difference of the leading physical anomalous dimension 5!). Indeed we 
have for £ DY — > 0: 

: (3.28) 



e DY tanh- 1 (Vl-4^ E 



Cdy ln(£ r 



Reporting the result (3.26) into (3.13), and performing the r DY — > expansion (r DY 



l-r 



under the integral (which is now legitimate, since the singularities of the integrand 



are all integrable) , one thus obtains the expansion of the physical evolution kernel: 

dp 2 



K dy (t,Q 2 )\ 



large fio 



C F 



{7- J 

' DY JO 



P 



pv(p 2 ) [MM-Mo) 



+ 



+ ^DY 



d^ 

p 2 

fOO 



Pv(p 2 ) (M(.M-M(P)) 



^Pv(p 2 ) (M(M - M(0)) + 0{rl Y )} 



(3.29) 



00 dp 2 



+ 5(1-t)C f I ^af"V)V t (e) , 
Jo A* 

where I performed integration by parts in the virtual contribution, in order to obtain at 
least a finite integrand, Vt(0) being infinite (the integral itself is infrared divergent since 
V t (0) = -1). Comparing with (3.10) I deduce (£ DY = p 2 /W 2 Y ): 



Si(W DY ) I large ^ = 

and also (e = p 2 /Q 2 ): 

dln\(n-Q 2 ))\ 2 



Cf r^-pvip 2 ) (fsAM-fsM) 
Jo P v 7 



(3.30) 



dlnQ 2 



large (3 



/°°^a^'" A (/r)V,(r) 
Jo 



A* 



(3.31) 
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which checks Eq. (3.10) in the large-/?o limit, and in addition gives the (large-/3o) dispersive 
representations of the physical Sudakov anomalous dimensions S and Si, as well as the 
(large-/?o) dispersive representation of the time- like quark form factor [14]. Finally, the 
dispersive representation of the regularized virtual contribution in (3.11) is given by: 



Jo A* Jo f 1 v 7 

(3.32) 



dlnQ 2 

where the integral on the r.h.s. is convergent both in the infrared and the ultraviolet region. 



3.3 Dispersive representation of the physical Sudakov anomalous dimensions 
(finite (3$) 

The generalization of Eq. (3.30) at finite /3q is straightforward [6,13]: 

roo j, ,2 



S{Wl Y ) = c F r^Lp s (^) (f s {M-Fs(oj) 

Jo ^ 7 

s t (wi Y ) = c F r** PSi tf) fe(e DY )-^(o)) , 

Jo f 1 v 7 



(3.33) 



where, in full analogy with (2.15), 

,Mink / , ,2 



( 2, da fX^) M in w 2, r dm2 ( ^ 

and, similarly to (2.16), ps,Si{l^ 2 ) correspond to the timelike discontinuities of some "Eu- 
clidean" effective charges, originally defined for spacelike momenta: 

Euclfl 2\ _ f Mink/ 2\ 

fJ- 2 /k 2 f°° dp 2 2 

From these dispersive representations, the perturbative expansions of the 'soft' Sudakov 
effective charges a^ rmk (p 2 ) and a^ mk (p 2 ) can be derived in a straightforward way [6] order 
by order in the full non-Abelian theory, given the expansions of S{W 2 Y ) and Si(W 2 Y ). 
Furthermore, the generalization of (3.31) is 

whereas that of (3.32) is: 

din | (H-Q 2 )) I 2 / Q2 rf/x 2 2 /■« d M 2 r Mmt 2 2 2 x 

— — + h l^ s{fi) - CF Jo /Q ) + ij (3 . 37) 

+ ar nfc (^)(^(^ 2 /Q 2 )-l)] , 
where terms have been arranged properly to have both infrared and ultraviolet convergence. 
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4. Conclusions 



I presented a simple momentum space ansatz for the structure of a systematic expansion of 
the physical evolutions kernels around the Sudakov limit, focussing on DIS and the Drell- 
Yan process. Each order in the expansion is given by a peculiar 'jet' or 'soft' (depending 
on the process) physical Sudakov anomalous dimension. 

The ansatz has been derived from a large-/3o dispersive representation, which can be readily 
extended to the full non-Abelian theory at finite Clearly more justifications and checks 
are still needed. In this respect, the general OPE method introduced in [3, 4] may give 
the appropriate tool for a systematic derivation. Another possibility (which shall be the 
subject of a future investigation) is to perform a check of the ansatz by matching it with 
existing [24, 25] fixed order calculations; as a by-product one would get (if the ansatz 
turns out to be correct) the perturbative expansion of the relevant 'jet' and 'soft' physical 
anomalous dimensions which occur as coefficients in these expansions. 
I have shown that the connection of the ansatz with the dispersive representation requires 
a non-trivial, and process-dependent, choice of expansion parameter in the Sudakov limit. 
While the adequate parameter r in the DIS case is standard, and refers to the final state jet 
mass scale, the meaning of the corresponding parameter r DY in the Drell-Yan case remains 
to be understood: it does not correspond (as could be naively expected) to the total energy 
E DY radiated in the final state. 

It was found that, although the leading jet {J{W 2 )) and soft {S{W 2 Y )) momentum space 
physical anomalous dimensions do not [6] have renormalons (at least at large-/?o), renor- 
malons do start to appear in the subleading jet (Jq(W 2 )) and soft {Sq{W 2 y )) physical 
anomalous dimensions. 

Given the multiplicity of emerging physical Sudakov anomalous dimensions, it would be 
interesting to find out whether there exists any relationship between them. In particular, 
one might wonder whether any exact, or approximate, universality property holds among 
the various 'Sudakov effective charges' aj in j_ and a^'^.. Actually, from the observation that 
the sum J'(W 2 ) + Jq{W 2 ) is a total derivative at large-/3o (which follows from the finitness 
of +-7\7 (£) at £ = 0, see (2.22)), and making the assumption that this property 

still holds at finite fio, one can already readily deduce in the DIS case that aj- ink (/x ) and 
aj^ k (fi 2 ) have identical expansions up to NLO, and therefore [6] coincide with the cusp 
anomalous dimension up to this order. 

The transition to an 1/N expansion in moment space has been shown to be straightforward, 
and the connection between the momentum and moment space jet (or soft) scales has 
been clarified. From the expansions in moment space, one can easily obtain the moment 
space dispersive representations, following the method in [6], by substituting the dispersive 
representations of the physical jet and soft physical anomalous dimensions in the moment 
space expansion. 

The longitudinal structure function Fl(x, Q 2 ), which has not been addressed in this paper, 
needs a special treatment (to be reported elsewhere). Indeed, it appears the dispersive ap- 
proach does not work in a straightforward way in this case, essentially because at large-/?o 
the longitudinal physical evolution kernel Kl(x,Q 2 ) does not coincide any more with the 
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derivative dF L (x, Q 2 )/dQ 2 . This is due to the fact that the longitudinal coefficient function 
is 0(a s ), rather then 0{oP s ) as in the case of F2 (which could be replaced at large-/?o by 
its leading term 5(1 — z) on the right hand side of Eq. (2.1)). In particular, Kl(x,Q 2 ) 
is not suppressed for x — > 1, contrary to Fl(x,Q 2 ) itself. Preliminary investigation never- 
theless seems to indicate that an ansatz of the form of Eq. (2.11) might still be valid for 
Kl(x, Q 2 ) (with a different coefficient of the (5(1 — x) term), although the connection with 
the dispersive approach is lost (except eventually at large-/?o, where the leading Sudakov 
physical anomalous dimension Ji (W 2 ^ could be a total derivative). 

Finally, although I have studied here the simplest examples where only one scale ('jet' or 
'soft') is present, the ansatz can be readily extended to other cases [6] where both of these 
scales do occur simultaneously. 
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A. N — > 00 expansion in moment space (Deep Inelastic Scattering) 

The N — ► 00 expansion of the moment space physical evolution kernel Eq. (2.3) can 
be straightforwardly obtained by taking the moments of the momentum space expansion 
Eq. (2.11) (r = ±=£ W 2 = rQ 2 ): 

K(N,Q 2 )= f 1 dxx N - ll -j{W 2 ) + dln if ( ^ 2 2))2 + C dxx^J^W 2 ) + ... (A.l) 
Jo r amQ J 

where the argument W 2 of the 'jet' anomalous dimensions, being x-dependent, has to 
be integrated over. The leading term in this expansion is infrared divergent at x = 1 
(r = 0), but, as we shall see below, this divergence is regularized by the virtual contribution. 
Considering first the subleading terms, which are infrared finite, we have identically (i > 0): 

£ d " H - 1 ' , *W = j^F* (t^t)" +1 ^.(4) < a - 2 » 

where I made the change of variable t = Nr. Then, using the N ^ 00 expansion (with 
p=l) 

( 1 \ N+P 1 t 

{tTJJ =e ~ t+ N e_ * 2 (t " 2P) + ° {1/N2) ' (A - 3) 

one gets, letting N —> 00 with the moment space 'jet' scale W 2 = Q 2 /N fixed inside the 
integral on the r.h.s. of (A. 2): 

P 1 "1 (* OO 

dxx N ' 1 r i J l (W 2 ) = J^TTlJ o dte-H'Ji {tW 2 ) + 0(l/N 2+i ) . (A.4) 

One should note that the integral on the r.h.s. of (A.4) contains constant terms, as well as 
logarithms of N. Moreover, the integral being convergent, the constant terms determine 
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the logarithmically divergent terms, given the beta function coefficients: see the argument 
following Eq. 4 in [13]— which implies that the knowledge of the logarithmic terms at a given 
subleading order of the 1/N expansion also fixes the constant terms at the same order. In 
particular, the N — > oo expansion of the next to leading real emission term in Eq. (A.l) is: 

J dx x"" 1 Jo {W 2 ) = — J dt e"* Jo (t W 2 ) + 0{l/N 2 ) . (A.5) 

Considering now the leading, infrared divergent term in Eq. (A.l), the same procedure 
yields the analogue of Eq. (A. 2): 

N+l 



f 1 dxx^ i -j(w 2 )= r 

Jo r Jo 



dt 



1 



1 + 



N 



1 



<2 Z 



(A.6) 



and, taking the N — > oo limit with Q 2 /N fixed inside the integral on the r.h.s., the analogue 
of Eq. (A.4): 

f 1 dx x N ~ 1 - J (W 2 ) = !°° dt e~* -j(t W 2 ) + 0(l/N) . (A.7) 
Jo r Jo t \ J 

The infrared divergence in Eq. (A.7) can be regularized as in Eq. (2.12) by subtracting from 
the r.h.s. of (A.7) the IR divergent piece ^rJ{n 2 ), and merging it together with the 

virtual quark form factor contribution m Eq. (A.l) (where the IR divergences 

cancell, as we have seen in section (2.1)). Since ^rJip?) = f Q N dt j J (t^-j (where I 



set t = Nfi 2 /Q 2 ), one obtains the N — ► oo expansion of the (regularized) leading term in 
Eq. (A.l): 



Jo 



dx x 



N-l 



;J(«'- 2 ) 



+ nJ ^e-*-(t-2)^(tW- 2 ) + 0(l/iV 2 ) . 



(A.8) 



Combining (A.5) with (A.8) one thus gets the large N expansion of the first two leading 
real emission terms in Eq. (A.l): 



dx x 



N-l 



J {W 2 ) + J (W 2 ) 
- + 



dt^\j(tw 2 )-l dt\j(tW 2 ) 

^ £ ' dte~ l (jo(tW 2 ) + \{t-2)j(tW 2 



+ 0{l/N 2 



to which 5 one has to add the (N-independent) regularized virtual contribution 
If fj(» 2 ) 



(A.9) 

d\n(F{Q 2 )) 



dhiQ 2 



+ 



5 The first term on the r.h.s of (A.9) can afso be written, up to exponentially small corrections at large 
N, as f N <tt(e-« -l)\j(tW 2 y 
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I also note that taking the moments of the leading term in the first line of Eq. (2.7) 
one gets (setting t = N(l — x)): 



f^-'T^-^H^H) M4) (A10 

= J dte- 1 -J (tW 2 )+0(l/N) , 
which, compared to (A. 7), checks the equivalence of the leading terms in (2.7) and (2.12) 



B. iV — > oo expansion in moment space (Drell-Yan) 

The N — > oo expansion of the moment space physical evolution kernel Eq. (3.3) can simi- 
larly be straightforwardly obtained by taking the moments of the momentum space expan- 
sion Eq. (3.10) (r DY = 2 W 2 Y = r 2 DY Q 2 ): 



K DY (N, Q 2 ) = jf 1 dr t n -' ± S {Wl Y ) + dln|( ^ Q f ))|2 S(l-r) + £ dr r N ~ l S (W 2 Y ) +... 

(B.l) 

The leading term in this expansion is infrared divergent at r = 1 (r DY = 0), but, as we 
shall see below, this divergence is regularized by the virtual contribution. Considering first 
the subleading terms, which are infrared finite, we have identically (i > 0): 



ly^rUS, (1*2,) = j^jf*. «** < B - 2 > 

where I made the change of variable u = N r DY , which, remarkably, leads to an integral of 
a form quite similar to the one on the r.h.s. of (A. 2). Then, using the N ^ oo expansion 
(A. 3) one gets, letting N — ► oo with the moment space 'soft' scale Q 2 /N 2 = W 2 Y fixed 
inside the integral on the r.h.s. of (B.2): 

P 1 1 P OO 

J dr r* DY S t (W 2 Y ) = ^ J du e~ u u* (u 2 W^ Y ) + 0{l/N 2+l ) . (B.3) 
Thus the N — > oo expansion of the next to leading real emission term in Eq. (B.l) is: 

P 1 1 P OO 

J dr t n - 1 So {Wl Y ) = -] du e~ u S (u 2 W^ Y ) + 0(1/N 2 ) . (B.4) 

Considering now the leading, infrared divergent term in Eq. (B.l), the same procedure 
yields the analogue of Eq. (B.2): 
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and, taking the N — > oo limit with Q 2 /N 2 fixed inside the integral on the r.h.s. of (B.5), 
the analogue of Eq. (B.3): 



f 1 dr r N - x — S (W£ Y ) =2 r du e~ u - S (u 2 C) + 0(1/N) . (B. 

JO r DY JO U \ / 



6) 



The infrared divergence in Eq. (B.6) can be regularized as in Eq. (3.11) by subtracting from 
the r.h.s. of (B.6) the IR divergent piece f® ^-S(u 2 ), and merging it together with the 

virtual quark form factor contribution — ^n~g2 ~ 111 Eq. (B.l) (where the IR divergences 



u 



N 2 



(where 



cancell, as we have seen in section (3.1)). Since ^i r S(fi 2 ) = 2 f Q N du^S 

I set u 2 = N 2 a 2 /Q 2 ), one obtains the N — > oo expansion of the (regularized) leading term 
in Eq. (B.l): 



/' 

Jo 



drr 



N-l 



' DY 



due~ u -S u 2 W, 



du-S[u z Wi 



f 

Jo 

1 du e"« \{u-2)S (u 2 W 2 Y ) + 0(l/N 2 



(B.7) 



Combining (B.4) with (B.7) one thus gets the large N expansion of the first two leading 
real emission terms in Eq. (B.l): 



Jo 



drr 



N-l 



s {wl Y ) 



+ So {W 2 Y ) 



due~ u -S (u 2 W 2 
u 



N 



1 

U 



du-S[u 2 W 2 



+ 



1 du e~ u (So (u 2 W 2 Y ) + \ (u - 2) S (u 2 W 2 Y ) 
+ 0(1/N 2 ) , 

(B.8) 



dln|(^(-Q2))| 2 



to which 6 one has to add the (N-independent) regularized virtual contribution 

if . 

I also note that taking the moments of the leading term in the first line of Eq. (3.7) 
one gets (setting u = N(l — r)): 



/■oo 

=2 / due~ u 
Jo 



(B.9) 



-S(u 2 Wi Y )+0(l/N) , 
u \ / 

which, compared to (B.6), checks the equivalence of the leading terms in (3.7) and (3.11). 



6 The first term on the r.h.s of (B.8) can also be written, up to exponentially small corrections at large 
N, as 2 f Q N du (e~ u - 1) ± S (u 2 Wl Y ) . 



- 20 - 



References 

G. Sterman, Nucl. Phys. B281 (1987) 310. 
S. Catani and L. Trcntadue, Nucl. Phys. B327 (1989) 323. 

R. Akhoury, M. G. Sotiropoulos and G. Sterman, Phys. Rev. Lett. 81 (1998) 3819 
[arXiv:hep-ph/9807330]. 

M. G. Sotiropoulos, R. Akhoury and G. Sterman, arXiv:hep-ph/9903442. 
R. Akhoury and M. G. Sotiropoulos, arXiv:hep-ph/0304131. 

E. Gardi and G. Grunberg, "A dispersive approach to Sudakov resummation," (to appear in 
Nucl. Phys. B) arXiv:0709.2877 [hcp-ph]. 

G. Grunberg, Phys. Rev. D29 (1984) 2315. 

Y. L. Dokshitzcr, G. Marchcsini and B. R. Webber, Nucl. Phys. B469 (1996) 93 
[hcp-ph/9512336]. 

S. Catani, Z. Phys. C75 (1997) 665 [hcp-ph/9609263]. 

W. L. van Neerven and A. Vogt, Nucl. Phys. B568 (2000) 263 [hcp-ph/9907472]. 
G. Grunberg, "Infrared finite coupling in Sudakov resummation," [hcp-ph/0601140]. 
G. Grunberg, Phys. Rev. D73, 091901 (2006) [hep-ph/0603135]. 
G. Grunberg, Phys. Rev. D74 (2006) 111901 [hep-ph/0609309]. 
S. Friot and G. Grunberg, JHEP 0709 (2007) 002, arXiv: 0706. 1206 [hcp-ph]. 
G. Grunberg, AIP Con}. Proc. 892 (2007) 268 [hep-ph/0610310]. 
S. Forte and G. Ridolfi, Nucl. Phys. B650 (2003) 229 [hcp-ph/0209154]. 
E. Gardi and R. G. Roberts, Nucl. Phys. B653 (2003) 227 [hep-ph/0210429]. 
A. Sen, Phys. Rev. D24 (1981) 3281; 
A. H. Mueller, Phys. Rev. D20 (1979) 2037; 
J. C. Collins, Phys. Rev. D22 (1980) 1478; 
L. Magnea and G. Sterman, Phys. Rev. D42 (1990) 4222. 
G. Curci and M. Greco, Phys. Lett. B 92 (1980) 175. 

D. Amati, A. Bassetto, M. Ciafaloni, G. Marchcsini and G. Veneziano, Nucl. Phys. B 173 
(1980) 429. 

S. Moch, J. A. M. Vermaseren and A. Vogt, Nucl. Phys. B688 (2004) 101 [hcp-ph/0403192]. 
J. A. M. Vermaseren, A. Vogt and S. Moch, Nucl. Phys. B724 (2005) 3 [hep-ph/0504242]. 
M. Dasgupta and B. R. Webber, Phys. Lett. B 382 (1996) 273 [arXiv:hcp-ph/9604388]. 



- 21 - 



